Why does air passage over forest yield more rain? 
Alternative interpretations of Spracklen et al. 2012 
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Abstract 

Spracklen et al. recently presented a pan-tropical study of rainfall and land-cover 
that showed that satellite-derived rainfall measures were positively correlated with the 
degree to which model-derived air trajectories had been exposed to forest cover. This 
result confirms the influence of vegetation on regional rainfall patterns suggested in pre- 
vious studies. However, we find that the conclusion of Spracklen et al. - that differences 
in rainfall reflect air moisture content resulting from evapotranspiration - appears un- 
dermined by methodological inconsistencies. We discuss some alternative explanations 
that require investigation. These include the distinct role of forest evapotranspiration 
in creating low pressure systems that draw moisture from the oceans to the continental 
hinterland. This alternative physical process is consistent with the empirical findings of 
Spracklen et al. but underlines a greater potential danger of forest loss than is suggested 
by their analyses of moisture recycling. 
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1 Introduction 



In their recent pan-tropical study ISpracklen et al.l |2012| examined how air exposure to for 



est cover influences subsequent rainfall from air moving over the tropical land surface. The 
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struction of the air trajectories from the observed wind fields and the use of daily rain- 
fall statistics - allowed them to analyze the apparent influence of the forest on rainfall 
on a short-term time scale. Such an approach apparently offers detailed insights into the 
vegetation-rainfall relationship. The quantitative concepts that underlie the interpretation 
of the obtained evidence justify a comprehensive a nalysis. 



In this note we revisit the arguments that led ISpracklen et al.l |2012| to conclude that 
the additional rain reflects a higher air moisture content resulting from forest evapotran- 
spiration. In Section 2 we show that the quantitative conclusion about the importance of 
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evapotranspiration varies with respect to the chosen time length of the air trajectories and 
the period of the recorded rainfalL In Section 3 we show that Spracklen et aL's estimate 
of post-deforestation reduction in Amazonian rainfall is similarly a function of the selected 
duration of the air trajectories. We argue that, dependent as they are on arbitrary choices, 
these particular conclusions do not offer quantitative informatio n about the actual atmo- 



spheric processes. The empirical findings of ISpracklen et al.l 2012| of a comparatively higher 



rainfall produced by forest compared to non-forest air are less affected by the choice of the 
key quantitative parameters on which their analysis was based. In the concluding section 
we discuss the physical variables that determine local rainfall. On this basis we suggest an 
alternative explanation to the observed patterns that can be tested with additional data. 

2 Additional rainfall and evapotranspiration 



Spracklen et al. |2012| used the Tropical Rainfall Measurement Mission (TRMM) data in the 
form of daily rainfall corresponding to 1° x 1° degree cells. They matched each local rainfall 
measurement to a backward air trajectory that described the air motion during the ten days 
pre ceding the measureme nt, Fig. [T] 



Spracklen et al.l |2012| quantified air exposure to forest by integrating the leaf area index 
(LAI) of the vegetation traversed over the 10-day air trajectory. This magnitude, denoted 
ELAI, is measured in time units (days), as LAI itself is dimensionless: 

L 



SLAI = LAI X ti, tL = - <tT = 10 d. 

u 



Here LAI is the mean leaf area index on the traversed territory, is time spent on the 
terrestrial part of the trajectory by the air moving with velocity m, L is length of the terrestrial 
part of the trajectory, tx is the total time length of the trajectory (including its oceanic part 
not shown in Fig. [1]). 

For each of the four tropical regions they studied, Spracklen et al.l 2012| divided all the 



10 days' trajectories into deciles of their ELAI. They found that rainfall produced in the end 
of trajectories with maximum ELAI ("forest air", top decile of SLAI) is several times higher 
than rainfall produced by the air with minimum ELAI (" non- forest " air, bottom decile of 
SLAI). For each trajectory they also calculated SET = ET x - the total amount of 
moisture acquired by model-derived evapotranspiration ET into the air as it moves from the 
beginning of the trajectory to the point of observation. Then they calculated the mean differ- 
ence AEET between the forest and non-forest air trajectories and compared AEET with the 
mean difference ARain in daily rainfall produced at the end of the trajectorie s, Fig. 1. Noting 



that AEET / ARain > 1 (see Table 2 in their Supplementary Information) ISpracklen et al 



2012| concluded that forest- derived evapotranspiration "more than accounts" for the addi- 
tional rain. 

But AEET (kg m~^) and the difference in daily rainfall ARain (kg m~^ d~^) are mag- 
nitudes of different dimensions (i.e. units). The value of AEET grows with increasing time 
length II of the terrestrial part of the air trajectory, while ARain is calculated for an arbi- 
trary chosen unit of time t/j = 1 d. Fig. [1] Being based on measures made over arbitrary 
periods, the ratio EET/ ARain does not carry information about atmospheric processes. 

3 Post-deforestation precipitation reduction 



Spracklen et~aD |2012] further used their findings to estimate how precipitation in the Ama- 
zon could be affected by large-scale deforestation. They calculated the dependence of rainfall 
on 10-davs' ELAI Fig, [lb . 



Spracklen et al.l 2012| calculated the likely LAI distribution in a deforested Amazon. As- 



suming that the air circulation does not change upon deforestation, they re-calculated ELAI 



2 



for all the trajectories. Then they recalculated the daily rainfall based on the established 
dependencies of local daily rainfall on SLAI. For example, if a trajectory was characterised 
by SLAI = 10 d but fell to SLA I = 3 d because of deforestation, then, according to the 

the rainfall at the end of the trajectory would reduce 
^Fig. [TJo. Having applied this procedure to the entire 
2012| concluded that deforestation would reduce rainfall 



analysis of iSpracklen et al.l |2012 



from about 9 mm d ^ to 4 mm d 



Amazon basin, ISpracklen et al. 



by 12% during the wet season and by 21% d uring the dry s e ason. They indicated that 
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this reflects the loss of precipitation recycling [Spracklen et al 
http:/ /www. nature.com/nature/journal/v489/n7415/full/naturel 1390.html] 

However, these estimates vary with respect to the chosen duration of the air trajectory 
tx (and of its associated terrestrial portion ti). Essentially, the choice of the trajectory 
duration determines the spatial scale at which remote deforestation influences the local 
rainfall observations. 

For illustration consider the forest trajectory in Fig. [T^. The air covers the distance from 
the coast to the point of observation in tj^ = 4 d and the forest has a mean LAI of LAI = 3, 
such that SLAI = 12 d for this particular trajectory. If half of the forest closest to the 
coast is now destroyed and turned to a desert with LAI = 0, then SLAI of our 4 days ' land 
trajec tory halves to SLAI = 6 d, see ([1]) . The dependence established by ISpracklen et al. 



[2012[, Fig. [TJo, predicts a significant post- deforest at ion reduction at the point of observation 
O because of the decline in SLAI. 

Now consider air trajectories with only t^, = 2 d. In such a case, deforestation of the 
coastal forest leaves unchanged the SLAI of the 2 days' air trajectory arriving at point O. 
Over the two preceding days the air will be moving over the continental interior forest and 
is not affected by the coastal deforestation. As SLAI has not changed, one would conclude 
that there will be no precipitation reduction at point O upon the same amount of regional 
deforestation. Thus, the estimated impact of regional deforestation on local precipitation 
depends on the duration chosen by the investigators for the considered air trajectories. 
However, quantitative estimates based on arbitrary choices do not have a predicatve power. 



4 Discussion: What determines rainfall? 

To understand how forest cover influences rainfall we must understand what physical pa- 
rameters determine rainfall. Rainfall occurs when moist air ascends and cools. From simple 
mass balance considerations precipitation rate P can be expressed as 

P = wq{l-^,hs), (2) 

where w, q and 7^ are, respectively, the upward air velocity, absolute humidity and water 
vapor mixing ratio at the level where condensation commences and 7c is mixing ratio at a 
height where it stops. In particular, if w; < (the descending air motion), any large-scale 
precipitation would be absent. 

Ignoring turbulent admixture of moisture into the air, moisture content in the point of 
observation go is equal to go = ge + SET — SP, Fig{T] It depends not only on the cumulative 
evapotranspiration along the trajectory SET = ET x t^, but also on the initial moisture 
content gs at the beginning of the terrestrial part of the trajectory and the cumulative 
terrestrial rainfall SP = P x ti th at depletes the air mo isture content. Remarkably, this 



rainfall was completely neglected by ISpracklen et al.l |2012| in their analyses, along with the 



ev aporation and precipita tion processes that occurred on the oceanic part of each trajectory. 



Spracklen et al.l |2012| noted that moisture content diminishes less in the forest than 
in the non-forest air from the beginning of the land part of the trajectory: (gfo — qm) > 
{QuO — qv&)- However, three additional requirements must be met in order to demonstrate 
that the additional rain is explained by the air's higher moisture content due to forest 
evapotranspiration. First, the intensity of convection (described by vertical velocity and the 
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completeness of condensation in the column) must be equal at the point of observation for 
the forest and non-forest air. Second, go must be higher in the forest versus non-forest air 
(O'fO > Qno)- Third, SET (and not or SP) must solely determine this difference. 

(Such an analysis would certainly be difficult. One set of complications are due to the 
recognised deficiencies of models in adequately representing the key processes of atmospheric 
moisture transport. Indeed, there is a signicant mismatch between the mode lled atmo s pheri c 



moisture convergence and the runoff observed in at least some studies [e.g., iMarengd . |2005|. 
It is a pertinent concern here that the moisture convergence for the Amazon basin derived 
from global circulation models is only half the figure actually observed. This implies that 
modeled parameters must be treated with caution and are likely to suffer major distortions. 
For example, while TRRM rainfall, wind directions and LAI are empirically observable vari- 
ables, evapotranspiration ET is model-derived. If the models underestimate the real moisture 
convergence but satisfactorily reproduce rainfall, this would mean that evapotranspiration 
is systematically overestimated.) 

Noting that any of parameters in ([2]) might account for the rainfall differences between 
the forest and non-forest air, we can propose a n alternative, and fully consistent, expla- 



nation for the empirical patterns established bv ISpracklen et al.l |2012| . Recently we have 



proposed that natural forest cover can lead to low atmospheric pressure. The mechanism 
derives from the subtle influences of atmospheric moisture, evaporation and condensation 
- in brief the areas with the highest evaporation drive upwelling and condensation, that 
induces low pressure and draws in most air from else where leading to a net atmosph eric 



moisture inflow to t he continent from the o cean (see [Makarieva and Gorshkovi . |2007| for 
the basic ideas and jMakarieva et al. . 2013b and references therein for a fuller account of 



the physical principles behind it). Evidence for these mechanisms has already come from a 
number of studies showing how rainfall over deforested areas tends to decrease exponentially 
with increasing distance from the ocean, while it stays more or less constant over forest s 



Makarieva and Gorshkovi . 120071 ISheil and Murdivarsol . I2009L iMakarieva et al.l . I2009L l2013a 



The empirical findings of ISpracklen et al.l 2012| are consistent with this interpretation 
and mechanism. Forest air will produce more rain because it is most often associated with 
low pressure systems (ascending air motion w > 0, high rainfall), while the air that has 
arrived from more arid regions should be more often associated with high pressure systems 
(descending air motion w < 0, low rainfall). This approach underlies a different role for forest 
evapotraspiration in driving rainfall patterns. The large-scale pressure gradients that drive 
the condensation-induced air motion are pr o portion al to the intensity of local condensation 



and, hence, precipitation [Makarieva et al.l . |2013b|. In the stationary case P = ET -|- C 



where C is the net amount of atmospheric moisture imported to the region. Rather then 
merely influencing the moisture content in the air that is passing over the forest, the process 
of evapotranspiration can impact regional atmospheric dynamics by enhancing rainfall and 
thus modifying the large-scale pressure gradients. This, in turn, enhances and stabilizes 
precipitation in a positive feedback loop. Conversely, decreased evapotranspiration would 
contribute to slowing down the regional atmospheric dynamics and making it less reliable. 
Another spectacular mechanism by which plants can influence condensation intensity and, 
hence, the pressure gradients is the biotically controlled emissi on of biogenic conden sation 



nuclei [e.g., Pohlker et al. . 2012 [. a process also mentioned by Spracklen et al. j2012 |. The 



key point is that deforestation can erode terrestrial low pressure zones on land and thus lead 
to a greater decline in rainfall than analyses of moisture recycling alone suggest. Interactions 
between regional pressure gradients and hydrological forest functioning await further studies. 

Acknowledgements. We are grateful to D. V. Spracklen, S. R. Arnold and C. M. Taylor 
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Figure 1: (a) Simplified scheme of analysis of ISpracklen et al.l |2012| . Two land trajec- 
tories, one over forest and one over non- forest, are illustrated to highlight key processes. 
Air with moisture content gfB or q^B at the beginning of the trajectory (B) moves for ti 
days over the forest or non-forested land, respectively, then comes to the point of obser- 
vation (O) where precipitation rates Pfo and PnO (kg m~^ d~^) are recorded. ETf and 
ETn, Pf and (kg m~^ d~^) are the mean evapotran spiration and p r ecipit ation rates 
along the forest and non-forest trajectories, respectively. ISpracklen et al.l 2012| calculated 
ASET = (ETf — ETn) xti and compared it with the difference in daily precipitation at the 
point of observation, ARain = (Pfo — Pno) x tR, tn = 1 d. Note that air moving at about 4 
m s~^ spends about 7 hr (much less than tfi) in the point of observation (1 x 1° cell, approx. 
100 X 100 km2). 

(b) The dependence of the observed wet season rainfall on SLAI established by 
Spracklen et al. j2012 | for Minas Gerais, Brazil. 
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